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Abstract: We consider the full Boltzmann equations for standard and soft leptogene- 
sis, instead of the usual integrated Boltzmann equations which assume kinetic equilibrium 
for all species. Decays and inverse decays may be inefficient for thermalising the heavy- 
(s)neutrino distribution function, leading to significant deviations from kinetic equilibrium. 
We analyse the impact of using the full kinetic equations in the case of a previously gen- 
erated lepton asymmetry, and find that the washout of this initial asymmetry due to the 
interactions of the right-handed neutrino is larger than when calculated via the integrated 
equations. We also solve the full Boltzmann equations for soft leptogenesis, where the lep- 
ton asymmetry induced by the soft SUSY-breaking terms in sneutrino decays is a purely 
thermal effect, since at T = the asymmetry in leptons cancels the one in sleptons. In 
this case, we obtain that in the weak washout regime {K < 1) the final lepton asymmetry 
can change up to a factor four with respect to previous estimates. 
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1. Introduction 

The discovery of neutrino oscillations makes leptogenesis a very attractive solution to the 
baryon asymmetry problem [1]. It is usually assumed that the tiny neutrino masses are 
generated via the (type I) seesaw mechanism [2] and thus the new singlet neutral leptons 
with heavy (lepton number violating) Majorana masses can produce dynamically a lepton 
asymmetry through out of equilibrium decays. Eventually, this lepton asymmetry is par- 
tially converted into a baryon asymmetry due to fast B ~ L violating sphaleron processes. 

Most studies of leptogenesis use the integrated Boltzmann equations to follow the evo- 
lution of the heavy particle number density and the lepton asymmetry. This approach 
assumes Maxwell-Boltzmann statistics, as well as kinetic equilibrium for all particles, in- 
cluding the heavy species. This assumption is normally justified in freeze-out calculations, 
where elastic scattering is assumed to be much faster than inelastic reactions. However, in 
the present context, kinetic equilibrium for the heavy species would have to be maintained 
basically by the decays and inverse decays alone, and it is not obvious that the integrated 
Boltzmann equation is always a good approximation. In general, 1 -^-^ 2 processes are more 
inefficient for thermalization compared to 2 2 processes, and in some parameter ranges 
there can be large deviations from kinetic equilibrium. 

In [3] , the impact of this difference on the lepton asymmetry produced during leptoge- 
nesis was studied, and it was found that in the strong washout regime the final asymmetry 
is changed by 15 — 30% when the full Boltzmann equations are used. In this work we 
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extend the study to two different scenarios not considered previously, in which the effects 
can be sizeable: 

(i) Preexisting lepton asymmetry. A lepton asymmetry that has been previously gen- 
erated, for instance by the next-to-lightest right handed neutrino, A?^2, tends to be washed 
out by the interactions of the lightest one, Ni. If M/Vi <^ 10^ GeV, the lepton asymmetry 
generated in its decay can be neglected and given that the A^i are thermally produced at 
temperatures close to Mi, this washout is exponential and therefore a change of order 20% 
may be important. We have used the full Boltzmann equations to calculate the evolution 
of the lepton asymmetry, created during N2 decay, at lower temperatures, T ~ Mi <C M2. 
We have also generalize them to the flavoured leptogenesis case. 

(ii) Soft leptogenesis. Since right-handed neutrino masses and therefore leptogenesis 
are usually associated to a very high energy scale, a supersymmetric scenario is desirable in 
order to stabilize the hierarchy between the leptogenesis scale and the clectroweak one. It 
has been shown [4-6] that supersymmetry-breaking terms can play an important role in the 
generation of a lepton asymmetry in sneutrino decays: they remove the mass degeneracy 
between the two real sneutrino states of a single neutrino generation, and also provide 
new sources of lepton number and CP violation. As a consequence, the mixing between 
the two sneutrino states generates a CP asymmetry in the decay, which can be sizable 
because of the resonant effect [7] of the two nearly-degenerate states. This scenario has 
been termed "soft leptogenesis", since the soft terms and not flavour physics provide the 
necessary mass splitting and CP-violating phase. It has also been studied in the minimal 
supersymmetric triplet seesaw model [8,9] and in the inverse seesaw scenario [10]. An 
important difference with respect to the standard leptogenesis mechanism, is that the 
lepton asymmetry produced through soft leptogenesis is a pure thermal effect, because at 
T = the asymmetry in Icptons cancels the one in slcptons. Only at finite temperature 
the difference between the fermion and boson statistics leads to a non-vanishing lepton 
and CP asymmetry, so in this case the use of the full Boltzmann equations is mandatory. 
Several approximations have been used in the literature [6,9], and we will compare our 
exact results with them. 

This paper is organized as follows. In section 2, we derive the full Boltzmann equations 
for standard leptogenesis and we study the washout of a previously generated lepton asym- 
metry. Section 3 is devoted to soft leptogenesis. In section 4, we present our conclusions, 
and more technical details concerning the full Boltzmann equations are described in the 
appendices. 

2. Standard leptogenesis 

In this section, we review the full Boltzmann equations relevant for leptogenesis, in a 
simplified model which includes only the heavy right-handed neutrino (iVj) decays, inverse 
decays, and resonant scattering. The off-shell 2 2 scattering processes mediated by Ni 
have only small effects for T < 10^^ GeV and can be neglected in first approximation [11,12]. 
We then investigate how the use of the full evolution equations affects the final lepton 
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na = ga I TTTZv? W ' (2-3) 



number asymmetry. We do not include thermal corrections [13] and we consider only 
initial zero abundance of the heavy neutrinos. 

The CP asymmetry in the decay of the right-handed neutrino Ni is: 

\A{N, ^ LH)\'' -\A{N, ^ LH)\^ ^ \A{Ni ^ LH)\^ - \A{Ni ^ Llif 
\A{Ni^LH)\^ + \A{Ni^LHf 

where we implicitly sum over all flavours, since, at this point, we work in the single flavour 
approximation. Its decay width is 

* a 

where l^j are the Yukawa couplings of the heavy neutrinos. 

We will denote fa the phase-space density of a particle species a, so its number density 
is given by 

where ga is the number of internal degrees of freedom. In order to eliminate the depen- 
dence in the expansion of the Universe, as usual, we write the equations in terms of the 
abundances, Ya = na/s, being s the entropy density, 

s = ^g.T^ (2.4) 

with the number of relativistic degrees of freedom at temperature T, so that = 106.75 
in the Standard Model (SM) and 3* = 228.75 in the Minimal Supersymmetric Standard 
Model (MSSM), which we consider in Sec. 3. 

We will study the time evolution of the right-handed neutrino distribution fjsf^, and 
the lepton asymmetry distribution fc = fL — fL- to the fast gauge interactions, to a 
good approximation the Higgs field and the leptons are in kinetic equilibrium. Moreover, 
the Higgs number asymmetry is not conserved due to the large top Yukawa coupling, so 
we neglect the Higgs chemical potential*. Then, we consider the following distributions: 

Z-' = (e^H/T (2.5) 

where we have introduced a chemical potential for the leptons, fi. Following [3], we use 
that /x/T <C 1 and we approximate: 

f^= ^e^./T + ,y f + <^iiffh (2.7) 
/^ + /i:^2/2^ + 0((M/r)2). (2.8) 



*The effect of keeping the Higgs number asymmetry has been studied in [14], and could lead to a 
reduction of the final baryon asymmetry of C(l). 
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where , 

/r = (e^-/^ + 1)-^ . (2.9) 
Then, the lepton asymmetry Yc = nc/ s \s given by 

Yc = n'^ + 0{{^f). (2.10) 

Note that in this section we have neglected the thermal masses mH{T), mL{T), therefore 
Eh,l = \Ph,l\ = Ph,l- a thorough study of thermal leptogenesis can be found in [13], where 
it has been shown that when thermal masses are taken into account, at sufficiently high 
temperature, the Higgs becomes heavier than iVj and the decay Ni HL is kinematically 
forbidden. At higher temperatures, the Higgs becomes heavy enough to allow the decay 
H Nj^L, where also a CP asymmetry en is produced. However this asymmetry turns 
out to have a negligible effect on the final results for leptogenesis, so we do not consider it 
here. 

Using the formalism of appendix A, we can write the full Boltzmann equations for the 
heavy neutrino and the lepton asymmetry. We assume that the later is small, so we work 
to first order in e and Yc- It is convenient to use the dimensionless variables Zj = Mi/T, 
Ei = EnJT and j/a = Pa/T (a = Ni,H,L) to eliminate the dependence on the expansion 
rate of the Universe, so we obtain ^: 

dfN, Kizf 



i \ dyn [/f /!^(1 - fNj - fN,{l - mil + r„^)] , (2.11) 



dzi ViEi 
dfc Kizf 



-Ifcifn' + fN,)] , (2-12) 

where we have defined the decay parameter Ki = Ti/H{T = Mi) , which controls whether 
or not Ni decays are in equilibrium. 

Equations (2.11) and (2.12) can be integrated numerically for given values of Ki and Cj. 
Note that within our first order calculation, in the evolution equation for the heavy neutrino 
distribution Eq. (2.11) only the equilibrium Higgs and lepton distributions appear, so it 
can be solved independently of the lepton asymmetry. Using the approximate expression 
Eq. (2.7) for fc, the integration of Eq. (2.12) over the dimensionless lepton momentum, 
UL, leads to the following evolution equation for the chemical potential, 

Td^i=^i ^y^|-r(rT^(^^+^^^^+ 

+e^ ifN^ - f2) [(1 - + /f ) - film } . (2-13) 



^Our equation for the evolution of fc is slightly different from the one in [3], due to the fact that they 
add the term 2ei(l — (in our notation) because of the resonant part of the LH ^ LH scattering. 
However, this resonant contribution involves f^, instead of /jv^ , as we show in the appendix, so we have an 
extra term '^^if'jff'j^ifNi — f^ )- have checked that this difference is numerically very small and does 
not change the main results. 
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Recall that /i is related to the lepton asymmetry Yc by Eq. (2.10). 

It is straightforward to verify that if Maxwell-Boltzmann statistics and kinetic equi- 
librium for all species are assumed, so that the phase space distributions are 

/iV. = Se-^S /^^ = e-^«, = e-^-±^/^ (2.14) 

the above equations can be easily integrated and one recovers the usual integrated Boltz- 
mann equations for leptogenesis: 



= -KMYN,-Y;.''y^:^, (2.15) 



^ = 6, KMYn, - Y^j'^^ - '1k,1CM)Yc , (2.16) 

azi l^2\Zi) 4 

where lC\{zi) and ]C2{zi) are the modified Bessel functions of the second kind of order 1 
and 2. 

We have solved Eqs. (2.11) and (2.13) in the case that the lepton asymmetry is gener- 
ated by Ni, the lightest right-handed neutrino. The results for the N\ abundance normal- 
ized to the equilibrium one Yn/Y^ are shown in Fig. II, left panel (blue solid line). As was 
already noticed in [3], at high temperatures the equilibration rate of the heavy neutrino is 
faster when the full Boltzmann equation is used, so the abundance is larger than the one 
obtained assuming kinetic equilibrium (pink dashed line). 

After conversion by sphaleron transitions, the resulting baryon asymmetry is related 
to the generated lepton asymmetry Yc by [15] 

Yb = ^Yc (2.17) 

In Fig. I we plot our results for the Baryon Asymmetry of the Universe (BAU) Yb 
(blue solid line) and the ones obtained using the integrated equations (pink dashed line) 

for Ki = 1, 10 and ei = 10~^. We find, in agreement with [3], that the difference between 
the two approaches is at most 20% for Ki > 1. For smaller values of Ki the difference can 
be larger, but in this limit our results are no longer valid because we have not included 
scattering which can enhance the neutrino production and change our results. 

Next, we focus on a different scenario: we assume that a sizeable lepton asymmetry 
has been produced initially (for instance, during N2 decay), and we study the washout of 
this asymmetry by Ni interactions, relevant at T ~ Mi <C M2, when the lepton asymmetry 
generated by A''i decays is negligible, which is always the case if M^Vj <C 10^ GeV [16,17]. 
Here, as discussed in [18], such lepton asymmetry may survive in two ways: (i) because A'^i 
interactions axe weak, so it only washes out a small part of the lepton asymmetry; (ii) if the 
A^i decays when only the r Yukawa interactions are in equilibrium, at 10^ GeV <T < 10^"^ 
GeV, generically a part of the preexisting lepton asymmetry is always protected from A^i 
washout due to flavour effects. We first study the effect of using the full Boltzmann equa- 
tions to follow the evolution of the preexisting lepton asymmetry, due to A^i interactions, 
in the unflavoured case. 
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Figure I: BAU obtained from the integration of the fuU Boltzmann equations, Eqs. (2.11) and 
(2.13), (blue sohd hne) and the approximate integrated equations, Eq. (2. 15), (2. 16) (pink dashed 
Unc) for Ki = l (left) and 10 (right), and ei = 10"^. 



We therefore solve eqs. (2.11) and (2.13) with the initial conditions at zi = Mi/T « 1 
of zero A'^i initial abundance and Yc = 10~^, for different values of the decay parameter 
Ki and ei = 0, i.e., we neglect the lepton asymmetry generated by Ni decays. In this ap- 
proximation, the Boltzmann equation for the lepton asymmetry can be solved analytically 
and we obtain: 



3Ki 



2 Pfi- 



dz^zf I + (2.18) 



We show our results in Fig. II. As we explained before, the A^i abundance (left panel) is 
independent of the lepton asymmetry (to the order we are considering), so it is the same 
as in the previous case. In the right panel we show the evolution of the baryon asymmetry, 
as a function of z = Mi/T. We obtain that the use of the full Boltzmann equations always 
decreases the final lepton asymmetry, as compared with the standard approximation. This 
can be understood from Eq. (2.18): in the usual Maxwell-Bolzmann approximation, only 
the distribution appears in the exponential washout factor. However in the full equation 
there is an extra term, /at., and thus the washout is stronger. We can see in this figure 
that for Ki < 1 the final asymmetry is only slightly reduced with respect to the integrated 
equations and this suppression increases with Ki. 

In Fig. Ill, we show the amount of initial lepton asymmetry that needs to be produced 
by in order to obtain the observed BAU after A'^i washout, as a function of the parameter 
Ki . The coloured region indicates the values of Y^" and Ki for which we obtain Yb within 
50% of the experimental value Yg = (8.75 it 0.23) x 10^^^ [19]. In the case considered 
here of hierarchical right-handed neutrinos, the asymmetry generated by is typically 
< iQ-8 Therefore, in typical scenarios, only for Ki < 3 we can expect a sufficient 
asymmetry to be generated. In fact, as shown in [17], we usually need Ki < 3 in 50(10)- 
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Figure II: Comparison between the integration of the full Boltzmann equations (Eqs. (2.11) 
and (2.13)) in N2 leptogenesis (blue solid line) and the approximate integrated equations 
(Eq. (2.15),Eq. (2.16)) (pink dashed line) for Ki = 0.1, 1, 10, and Yc = 10"^. In the plots on 
the left we show the heavy neutrino (iVi) abundance as a function of z, and in the ones on the right 
the surviving Baryon Asymmetry after A^i washout. 

inspired scenarios with flavour effects. 

Finally, we consider the addition of flavour effects in this scenario of an initial lepton 
asymmetry and its washout by A^i interactions. As shown in Appendix A, the Boltzmann 
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Figure III: The coloured region sliows the amount of initial lepton asymmetry that needs to be 
produced by N2 in order to obtain the right amount of BAU after iVi washout, as a function of the 
washout parameter Ki. 



equation for the A^i abundance is identical to the single flavour case, while now we have to 
consider different evolution equations for the asymmetries in the different flavours: 




[fH'fi'i^ - fN.) - fN.a - /r)(i + m ' (2-19) 

{K^ef (/^, - Q) [(1 - /f )(1 + /^^) - /^Vr] - 



-^fcMn' + /ivj} + 0{Y^.,Y^., X fc,) , (2.20) 

where K" = T(Ni — > Lo,H,LaH)/H{T = Mi). In this expression we can see that, up to 
subleading corrections in the Yukawa couplings (barring special cases where Y^-fc^ > /£„)> 
the evolution equations for the asymmetries in different flavours decouple and the situation 
is analogous to the single flavour case. The only difference in this equation is the presence 
of different CP asymmetries in different flavours and the fact that the washout is produced 
only by the inverse decays in the relevant flavour [16]. Therefore, we can extend the results 
obtained above in the single flavour approximation for instance to the three flavour case 
considering independently the asymmetries in the e, fi and r channels with independent 
washout Kf, and Kl. Figs. II and III are also valid for the asymmetries in the 
different flavours with the corresponding washout parameter. Although our results are 
model independent and can be used in a generic model, they basically agree with the 
results obtained in Ref [17] in an 50(10) scenario with integrated equations as they have 
always small washout in the r-flavour in the relevant parameter space. 
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3. Soft leptogenesis 

In this section, we briefly review the soft leptogenesis scenario, and we present the full 
Boltzmann equations relevant for soft leptogenesis in the context of the supersymmetric 
type I seesaw model. We solve them exactly and compare with the approximations used 
in the literature. 

For a hierarchical spectrum of right-handed neutrinos, successful leptogenesis requires 
generically quite heavy singlet neutrino masses [20], of order M > 2.4(0.4) x 10^ GeV for 
vanishing (thermal) initial neutrino densities [20,21], although flavour effects [11] and/or 
extended scenarios [22] may affect this limit ^ . The stability of the hierarchy between this 
new scale and the electroweak one is natural in low-energy supersymmetry, but in the 
supersymmetric seesaw scenario there is some conflict between the gravitino bound on the 
reheating temperature and the thermal production of right-handed neutrinos [24]. This 
is so because in a high temperature plasma, gravitinos are copiously produced, and their 
late decay could modify the light-nuclei abundances, contrary to observation. This sets an 
upper bound on the reheating temperature after inflation, Trh < 10^"^ GeV, which may 
be too low for the right-handed neutrinos to be thermally produced [24] . 

Once supersymmetry (SUSY) has been introduced, leptogenesis is induced also in 
singlet sneutrino decays. If supersymmetry is not broken, the order of magnitude of the 
asymmetry and the basic mechanism are the same as in the non-supcrsymmctric case. 
However, as shown in Refs. [4-6], super symmetry-breaking terms can play an important 
role in the lepton asymmetry generated in sneutrino decays because they induce effects 
which are essentially different from the neutrino ones. Soft supersymmetry-breaking terms 
involving the singlet sncutrinos remove the mass degeneracy between the two real sneutrino 
states of a single neutrino generation, and provide new sources of lepton number and CP 
violation. As a consequence, the mixing between the two sneutrino states generates a 
CP asymmetry in the decay, which can be sizable for a certain range of parameters. In 
particular, the asymmetry is large for a right-handed neutrino mass scale relatively low, 
in the range 10^ — 10^ GeV, below the reheating temperature limits, what solves the 
cosmological gravitino problem. This is the so called "soft leptogenesis" scenario, which 
we are going to consider now. 

The superpotential of the supersymmetric seesaw model contains the following relevant 
terms: 

W = ^MijNiNj + YaiNi (LaH) , (3.1) 

where the parentheses indicate SU(2) contractions and a,i = 1,2,3 are flavour indices. 
La,Ni,H are the chiral superflelds corresponding to the left-handed lepton doublets, the 

right-handed (RH) neutrinos and the up-Higgs doublet, respectively, and Y^i denote the 
neutrino Yukawa couplings (notice that we are using the convention a = L, j = R). The 

*This bound applies when the lepton asymmetry is generated in the decay of the lightest right-handed 

neutrino. The possibility to evade the bound producing the asymmetry from the second lightest right- 
handed neutrino has been considered in [23], and flavour effects have been analysed for this case in [16]. 
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soft supersymmetry-breaking terms involving the heavy sneutrinos are 

1 



C 



soft 



-mlN*N, 



AaiYaiNiiLaH) + -BijMijNiNj + h.c. 



(3.2) 



Contrary to the traditional leptogenesis scenario, where at least two generations of RH 
neutrinos are required to generate a CP asymmetry in neutrino/sneutrino decays, in this 
mechanism for leptogenesis, a single generation of heavy RH neutrinos is sufficient to 
generate a CP asymmetry in sneutrino decays. Therefore, from now on we consider a 
simplified one-generation model, which refers to the lightest of the three heavy sneutrinos, 
that we denote as 1. For simplicity we also assume proportionality of soft-trilinear terms, 
and drop the flavour index for the coefficent A. The sneutrino interaction Lagrangian is 
then: 



C 



(m^ + \M\^)N*N 



BMNN + h.c. 



YalN{Lah) + MYalN*{LaH) + AYalN{LaH) + h.c. 



(3.3) 



where h is the fermionic partner of the Higgs doublet H. Under these conditions, a physical 
CP-violating phase is still present in the neutrino sector. 



$ = arg{AB) 



(3.4) 



which we choose to assign to A. The right-handed neutrino has a mass M, while the 
sneutrino and antisneutrino states mix in the mass matrix, with mass eigenvectors 



iV+ = (e^*/2iV + e-**/2iV*) 
^*/2 jv - e-^*/2 TV*) 



and mass eigenvalues 



(3.5) 



(3.6) 



We define the fermionic and scalar CP asymmetries in the decay of each Ni (i = ±) as: 



\A{Ni^LH)\^-\AiNi^L^H^)\^ ^ \A{Ni ^ LH)\^ - \A{Ni ^ L^H^)\' 
\A{N,^LH)\^ + \AiNi^L^W)\^ 
\A{Ni^Lh)\^-\A{Ni^Lh)\-' \AiNi 



s\2 



■(3.7) 



"•^^ \A{Ni^Lh)\^ + \A{Ni^Lh)\^ 
where we have implicitly summed over flavours, and 

\Alf = 2 j2\yaiM\\ 

a 

M? 



Lh)\^ - \A{Ni ^ Lh)\' 



\A 



f\2 



(3.8) 



\A{f = 2j2\yaiM\'j^{l-XL-Xf,), 



(3.9) 
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with Mi = M+,M_ and 

(3.10) 

Notice that, in this section, we keep the thermal masses of the sneutrino decay products, 
since they break supersymmetry and contribute to obtain a non-vanishing CP asymmetry. 
They are given by [13]: 

mUT) = 2mUT) = (^^i + '-g^y + ^y,^) , (3.11) 

m|(r) = 2mi(r) = g^i + Ig'y) , (3-12) 

where 52 and gy are the gauge coupUngs and Yt is the top Yukawa. Then, the total 
sneutrino decay width is 

Ti = ,^^{X{1,x~^,xh)\AI\'' + \{1,xl,xh)\A{\^) , (3.13) 

where 

A(l, X, y) = ^/{l + x-y)^-Ax (3.14) 

The mixing between the sneutrino states can generate a sizable CP asymmetry in their 
decay, due to the resonant enhancement of the self-energy contribution. The CP asymmetry 
can be computed following the effective field-theory approach described in [25], which takes 
into account the CP violation due to mixing of nearly degenerate states by using resumed 
propagators for unstable (mass eigenstate) particles. Neglecting supersymmetry breaking 
in vertices and keeping only the lowest order contribution in the soft terms, it is given 
by [4,6]: 

4rB_W 

Here, we have neglected thermal corrections to the CP asymmetry from the loops, i.e., we 
have computed the imaginary part of the one-loop graphs using Cutkosky cutting rules 
at r = 0. These corrections are the same for scalar and fermionic decay channels, (only 
bosonic loops contribute to the wave-function renormalization common to both decays), so 
they are not expected to introduce significant changes to our results. 

In order to generate enough asymmetry, we need i? ~ F, thus the lepton number 
violating soft bilinear coupling B, responsible of the sneutrino mass splitting, has to be 
unconventionally small. Moreover, as one can see in the above equation, in soft leptogenesis 
induced by CP violation in mixing, an exact cancellation occurs between the asymmetry 
produced in the fermionic and bosonic channels at T = 0. Therefore, thermal effects play 
a fundamental role in this mechanism: final-state Fermi blocking and Bose stimulating 
factors, together with the effective masses of the particle excitations in the plasma, break 
supersymmetry and remove this degeneracy. 

Several comments are in order: 



^Neglecting soft SUSY-breaking corrections and thermal masses \A{f = \Af 
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The effects of flavour in soft Icptogenesis have been studied in [26], and we will not 
consider them here, since the single-flavour approximation is enough to illustrate our re- 
sults. 

It has been recently pointed out that for resonant scenarios the use of quantum BE may 
be relevant [27] . For standard resonant leptogenesis [7] , they induce a T dependence in 
the CP asymmetry which can enhance the prodiiccd baryon number. However, in Ref. [28] 
it has been shown that in soft leptogenesis, due to the thermal nature of the mechanism 
already at the classical level, the introduction of quantum effects does not lead to such 
enhancement and therefore it is enough to consider only the classical BE in this work. 

The authors of Ref. [5] identified new sources for soft leptogenesis, induced by CP 
violation in right-sneutrino decay and in the interference of mixing and decay. These 
contributions are relevant both because they can be sizable for natural values of the B 
parameter and because, unlike the CP violation in mixing, they do not require thermal 
effects, as they do not vanish at T = 0. However, this calculation has been recently 
revisited in Ref. [29], where it has been found that for all soft SUSY breaking sources of 
CP violation considered, at T = the exact cancellation between the asymetries produced 
in the fermionic and bosonic channels holds. Therefore it seems that the full Boltzmann 
equations are always required to calculate the final lepton asymmetry generated in soft 
leptogenesis. 

3.1 Full Bolzmann equations 

We assume that the sneutrinos are in a thermal bath with a thermalization time F"-^ 
shorter than the oscillation time, (AM)~^, therefore coherence is lost and we can write 
the Boltzmann equations for the mass eigenstates Eq. (3.5). As in the previous section, we 
work in a simplified scenario including only sneutrino decays, inverse decays and resonant 
scattering. Within this approximation, we can neglect RH neutrino interactions, since in 
soft leptogenesis only sneutrino interactions generate the lepton asymmetry. 

We assume that, because of the fast gauge interactions, the Higgs and higgsino fields 
are in thermal equilibrium and the leptons and sleptons are in kinetic equilibrium, so that 
the corresponding distributions are: 

•^"^^ exp[(Ei-M/)/r] + i ' •^^ = exp[(£;i + /x;)/r] + i ' ^^-^^^ 

^^ = exp[(Sj-M.)/r]-l ' ^^^ = exp[(Ej + /x,)/T]-l • ^^'^^^ 

We have introduced a chemical potential for the leptons, and sleptons, jig- We are 
interested in the evolution of the sneutrino density distributions /^^ and the fermionic and 
scalar asymmetries, /r = /l — /l, which is given by Eq. (2.7) to first order in /x/, and 

2 pEh/T .. .. 
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We approximate fL + f^c^ 2/^^ and fj^ + fj^^ ~ ^>^Z^' "^^6^6' 

In the limit niaiT) <C T, the chemical potentials are related to the corresponding asym- 
metries by ^ 

Yc = ,if^ + 0{{!^)'), (3.21) 



^£ = ^ + • (3.22) 



Note that the masses and widths of the two sneutrino states are equal as long as we 
neglect supersymmetry breaking effects, then ~ In ^ ^jv' ^^"^ '^^'^ write a unique 
BE for f^; thus, in the following, we do not write the subindex i. The total fermionic (or 
scalar) asymmetry is then twice the asymmetry generated by one of the two sneutrinos. 

The relevant Boltzmann equations are derived in appendix B, and can be written 
in terms of the dimensionless variables z = M/T, Ea = Ea/T and ya = Pa/T (a = 
iV, h, L, H, L) as: 



Z2 



Kf I dE, [fl^fl\l + /^) - /^(l - )(1 - )] + 

7 EI" 



dz 2yNEN 

+ Ks P dEn [/^Vf (1 + fri) - + f^)i^ + )] [ . (3-23) 

J E'^ J 



1 1 1 

dfr KsZ 



^fc{f^ + fl')\ + ^/,, (3.24) 



1.,. ,e,^ 1 



where 



+ 



+ oWn, - fn') \ + H-z^9^ (3.25) 



/f = (e^^/"^ - 1)-' , (3.26) 



and we have used the notation 

Kf = {I- XL- Xh)\{l,XL,Xh)K , Kg = X{l,XL,Xh)K , (3.27) 



^Taking into account the (s)lepton thermal masses mL{T) = 0.3T and mj^{T) = 0.4T, to first order in 
Hf,l^s we obtain 

/T~i2 rTi'2 

Yc = M/ 0.329 — , Y£ = /i, 0.541 — , 

which are the actual values used in our calculation, instead of the 1/3 and 2/3 of Eqs. (3.21) and (3.21), 
respectively. We see that the lepton thermal mass does not change significantly the lepton asymmetry, but 
the slepton thermal mass leads to a 20% reduction of the slepton asymmetry, proportional to 0.54 instead 
of 0.67. 
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with K = r^/H{T = M) defined in terms of the sneutrino decay width without thermal 
masses, 

r° = ^El^"^l^- (3-28) 

a 

Since the thermal masses are different for the final states hL and HL, so are the 
integration limits, given by: 

'"^ = i {Ei{l -XL + Xh) ± Vi A(l, XL,Xh)} , (3.29) 

f _ EL + 4zi} + XL-Xh)\{l,XL,Xh) 

'~ Il±^ + feA(l,x.,..) ' ^'-'^^ 

and analogously for E^'"^ and Ef, just replacing L, h with L, H in the above equations. 

The terms Sg and Sg in Eqs. (3.24) and (3.25) represent the fast gaugino interactions, 
defined in Eq. (B.ll). These interactions are in equilibrium and mediate processes that 
transform leptons into scalar leptons and viceversa [L + L L + L). Thus we shall impose 
that /U/ = 

As in the previous section, decays, inverse decays, and on-shell scattering processes 
must be considered in order to obtain the appropiate out-of-equilibrium condition (see 
appendix B). Using the approximated distributions (2.7) and (3.19) for fc and re- 
spectively, and integrating over the dimensionless (s)lepton momentum, y^, we obtain 
the following Boltzmann equations for the chemical potentials, defined in Eqs. (3.21) and 
(3.22): 



T dz 



— '^^y^j, '^^^Tp^^^^+^'^)- 
- ^ (/iv - /f) [(1 - - /D + /r/r] } - (M/ - ^^s)l, , (3.31) 



+ ^ (/iv - /f) [(1 + /?)(! + ) + fifn] } + (m/ - Ms)% • (3.32) 



where we have summed over the two sneutrino states and the integration limits are: 

1 
2 



E^'"' = l{ENil + XL-Xh)±yNXil,XL,Xh)} , (3.33) 



and analogously for E^'"^, just replacing L, h with L, H in the above equation. The jg, jg 
terms in Eqs. (3.31) and (3.32) represent gaugino mediated processes L + L -ir^ L + L. 
Schematically, we can write the equations for the chemical potentials as: 



dz 
dus 
dz 



eA + CfXf-{fif-HshgT , (3.34) 
eB + Dns + inf- l^s)lgT , (3.35) 
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while the lepton and slepton number asymmetries are related to the chemical potentials by 



Yc = afHf = asHs ■ (3.36) 

Therefore, the evolution equation for the total lepton asymmetry, Y^^ =Y£ + Yg, is 

dYr 

-j^ = eiafA + asB) + afCnf + asDns (3.37) 

Now, in order to take into account that at shorter time intervals the fast gaugino interac- 
tions make //^ = we replace this equality in the right-hand side of the above equation, 
obtaining 

= e {a, A + asB) + ^l£±^Ycr (3.38) 

We have solved numerically Eqs. (3.23) and (3.38). Following previous approaches [6], [26] 
we neglect the thermal masses in the evolution equation for the sneutrino distribution, 
Eq. (3.23), so the decay channels are always open. If one keeps the thermal masses, there 
is a range of T for which the decays N ^ hL, HL are kinematically forbidden [13] . However 
we have checked that for the relevant values of z, the effect of the thermal masses in the 
sneutrino distribution is negligible. 

It is important, though, to keep the thermal masses in the evolution equation for 
Yc^, both, because they contribute significantly to the asymmetry and because the Bose- 
Einstein distribution is divergent at low values of z for massless scalars, leading to an 
unphysical enhancement of the slepton number asymmetry. In this case, one should also 
consider the CP asymmetry produced in the decays H NiL and h NiL, allowed at 
higher temperatures. In [13], this contribution has been found to be negligible in standard 
leptogenesis, so we assume that this is also the case in soft leptogenesis. 

3.2 Approximated Boltzmann Equations 

The standard integrated equations, Eq. (2.15) and Eq. (2.16), can not be used in the 

case of soft leptogenesis, since the CP asymmetry produced in fermionic decays is exactly 
canceled by the one produced in the scalar channel if Maxwell-Boltzmann statistics is 
assumed. However one can find some approximate solutions in the literature which try to 
estimate the baryon asymmetry generated in soft leptogenesis. 

One possibility, which was used in [6] and [26], is to neglect all the thermal correc- 
tions except for the ones that are crucial to get a non vanishing CP asymmetry, which 
are evaluated in the approximation of decay at rest of the heavy sneutrinos. Then, the 
Boltzmann equations arc integrated in the standard way, i.e., assuming kinetic equilibrium 
and Maxwell-Boltzmann statistics for all the particles in the plasma, obtaining: 

^^ = -^.(r,->^')|W, ^ (3.39) 

= 2.(T)if.(y^-yi«)|M-:^c,(.)F.,, (3.40) 

where Ycj, = Yc + Y^, and we have summed over the two sneutrino states. 
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The effective, temperature-dependent CP asymmetry, e(T), includes the statistical 
factors and thermal corrections which are different for fermions and scalars, and it is 
defined as: 

e{T) = e '-^^ , (3.41) 

^ ^ CB + CF ^ ' 

where 

CB = X{1,xi,xh)[1 + fB{Ei)][l + fsiEH)] , (3.42) 

CF = {I - XL - Xh)X{l,XL,Xh)[l - fF{ELm - fF{Eh)] , (3.43) 

being /b and /f the Bose-Einstein and Fermi-Dirac distributions, respectively, and 

— z — z 

^L,H= 2^'^ + ^l,H- ^H,l) ' EL,h= ^il+XL,h- Xh,L) (3.44) 

We expect this approximation to be accurate in the strong washout regime, K ^ 1, 
since in this case the kinetic equilibrium for the heavy sneutrino is a good approximation 
and, moreover, the final lepton asymmetry is independent of the initial conditions, and 
fixed only by the late time z > 1 (low temperature) evolution, when the thermal motion 
of the heavy sneutrino can be neglected. We discuss the details in Sec. 3.3. 

We have also explored a different approach, based on Ref. [9]: we again assume kinetic 
equilibrium for all the particles in the plasma, as well as Maxwell-Boltzmann statistics for 
the heavy sneutrino, but we keep the phase space and statistical factors which are crucial 
in soft leptogenesis. Furthermore, we approximate the Fermi and Bose distributions by the 
Maxwell-Boltzmann one: 

/^(B) = ;^7^ - e-^/^ (3-45) 

since this is enough to have a non vanishing total CP asymmetry in sneutrino decay, and 
the neglected terms are further suppressed by extra 0{e~^^'^) factors. 

Taking into account the (s)lepton thermal masses, given in Eq. (3.12), the relation 
between the (s)lepton asymmetry Yc and the corresponding chemical potential is now 

= /x/^/C2(m^/T) + 0{{f^f) , (3.46) 

and the same expression holds for Y^, just changing /if ^ fig and L ^ L. 

Within these approximations, the Boltzmann equation for the heavy-sneutrino abun- 
dance is again the standard one, Eq. (3.39), while the evolution equations for the (s)lepton 
asimmetries become: 

dz 27r^ Jey^ 2 m|/C2(mj/r) L 

+ ^ifN- O [(1 + /?)(! + fH) + f^fS] • (3-48) 
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The integration limits E^'-^ are given in Eq. (3.33). Recall that fast gauge interactions 
imply that jif = jJLs- 

The integration over the (s)lepton energy can be performed analytically, and by writ- 
ting the result as a series expansion in the heavy-sneutrino momentum, j/jv, it is also 
possible to compute the integral over E'^r, order by order inyj^. The calculation is lengthy 
but straightforward, and it is described in appendix C. The integrated Boltzmann equation 
for the total lepton asymmetry Yc^ can be written as: 

^^ = 2eK{Y^- rp^^ - ^K:,{z)F,{z)Ycr ■ (3.49) 

Fi{z) is a complicated function of the thermal masses that can be found in appendix C. 
We present here the limit of vanishing thermal masses, for the purpose of comparison with 
Ref. [9]: 



^r(^)=4^ 3n+l22n-l,, ^n+l(3./2) 

n=0 

R 2z z 

~ -/Ci(3z/2) + y/C2(3z/2) + -/C3(3z/2) + . . . , (3.50) 

where the dots stand for negligible contributions. The approximation used in [9] seems to 
correspond to just keeping the first term in Eq. (3.50), although these authors obtain that 
the argument of the Bessel function K\ is \/2z instead of 3z/2. We find that to keep only 
the first term in F\{z) is not enough to get an accurate result. 
The function F2{z) is given by 



^ . ^ _ o^2 [A(l..>-L..r/,)]-(l - X, .ru) + [A(l. .r^. .r,,)]- 



It is easy to see that it reduces to the standard washout term in the limit of zero thermal 
masses, since /C2(m/T) ~ 2T^/m^ for m/T <^ 1. 

We expect that this approximation works better than the effective CP asymmetry 
e(T), since here the thermal motion of the sneutrino is fully accounted for. In fact, we have 
checked that our approximated equation (3.49) with the first three terms of Fi (z) gives the 
same results as the numerical integration, when one assumes kinetic equilibrium for the 
heavy sneutrino and Maxwell-Boltzmann statistics for all the particles. 

3.3 Results and discussion 

Our main results in the comparison between the three approaches to the Boltzmann equa- 
tions in soft leptogenesis are summarized in Fig. IV for different values of K. In the left 
panel, we plot the sneutrino number density normalized to the equilibrium density, from 
the full (solid blue) and integrated (dashed pink) Boltzmann equations. In the right panel, 
we show the baryon asymmetry obtained by integrating the full set of Boltzmann equa- 
tions, Eqs. (3.23) and (3.38), (solid blue), our approximated equations, Eq. (3.49) (dashed 
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Figure IV: Integration of the full Boltzmann equations in soft leptogenesis and the approximate 
equations, for K = 0.1, 1, 10, and e = 4.4 • 10~^. Left panel: sneutrino distribution as a function of 
z, obtained from the integration of Eq. (3.23) and the approximate Eq. (3.39). Right panel: baryon 
asymmetry obtained from the integration of Eqs. (3.31) and (3.32) and the generated in the two 
approximate equations Eq. (3.49) and Eq. (3.40). 



pink), and the effective e(T) approximation, Eq. (3.40), (dotted light blue). Recall that in 



-18- 




Figure V: Sncutrino momentum distribution assuming kinetic equilibrium (left) and solving the 
full Boltzmann equation (right), for different values of z, and K — 0.1, 1, 10. 



the MSSM, the final baryon asymmetry induced by sphaleron transitions is 

Yb = ^Yc, . (3.52) 
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We find that at high temperatures (small z), the equilibration rate of the sneutrino 
number density is higher when the full Boltzmann equations are used, similarly to what 
we have in Figure II in the case of the heavy neutrino. The main reason for this effect is 
that when z < 1, the low momentum states of the sneutrino are populated very efficiently, 
while the population of high momentum modes is small and comparable to the kinetic equi- 
librium distribution. This can be seen in Fig. V. where wc plot the sneutrino phase-space 
distribution, assuming kinetic equilibrium (left panel) and using the full Boltzmann 
equations (right panel), for different values of z. From these plots, it is clear that, even 
ioT z > 1 in the weak washout, the density of low momentum states is significantly larger 
than the density obtained assuming kinetic equilibrium. 

The behaviour of the baryon asymmetry, shown in the right panel of Fig. IV, results 
from the competition of different effects. In soft leptogenesis, due to the inclusion of the 
thermal masses, the CP asymmetry vanishes for z < 0.8, because both, fermionic and 
bosonic sneutrino decay channels are kinematically forbiden. For z < 1.2 the fermionic 
channel of sneutrino (inverse) decay creates an asymmetry, which tends to flip sign when 
the bosonic channel is open {z > 1.2), since this one dominates. On the other hand (as 
it happens also in the standard leptogenesis case), the asymmetry generated during the 
production has opposite sign (we call it "wrong-sign" asymmetry) to the one produced 
at later times, in N decay ("right-sign" asymmetry). In the different panels of Fig. IV 
one can see the competition of these two effects, depending on the strength of the Yukawa 
sneutrino interactions. 

In the case K = 0.1, we observe (in agreement with [6]) the flipping of the lepton 
asymmetry sign when the bosonic channel is open, both in the full (solid blue) and in 
the approximated calculations. Then, since the Yukawa interactions are weak and the 
decay occurs at late time, the generation of the right-sign asymmetry is not enough to 
overcome the wrong-sign one, created during the sneutrino production. However there is an 
important difference between assuming or not kinetic equilibrium for the heavy sneutrino. 
As can be seen in the first row of Fig. V, the sneutrino momentum distribution obtained 
with the full equations and for all values of z is always larger than the kinetic-equilibrium 
distribution. From Eqs. (3.24) and (3.25), we can sec that the source of the asymmetry is 
always proportional to — f^, where is the Bose-Einstein distribution. In Figure VI, 
we can see that up to 2; < 1, — ~ while for 2; > 1, — > for most values 
of the momenta. As a consecuence, the right-sign asymmetry for z > 1 is larger than in 
the kinetic-equilibrium case, partially compensating the wrong-sign asymmetry generated 
in the production although it is still not big enough to flip the sign again. Notice that 
this is not the usual washout, which is very small for K = 0.1. If one assumes kinetic 
equilibrium, the population of the momentum states is more uniform, and the right-sign 
asymmetry generated is smaller, so one gets a larger final result. However for ii' < 1 one 
has to be cautious, since scattering processes, that we have not included, can be relevant 
(see for instance [30]). 

For = 1, we see in the second row of Fig. IV that assuming kinetic-equilibrium 
distributions we find two sign-flips: one due to the opening of the bosonic decay channel 
and the second, at later time, because for larger K, the washout of the initial wrong- 
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Figure VI: Sncutrino momentum distribution normalized to tlic equilibrium one, the Bose-Einstein 
distribution. We show the ratio fjy/ for different values of z and A' = 0.1 (left), K — 1 (right). 

sign asymmetry is more efficient, and an asymmetry of the right-sign is finally created. 
However, using the full Boltzmann equation for the sneutrino, these two effects partially 
compensate each other. In fact, the bosonic channel is open at z > 1.2 and, as we see in 
Figure VI, — > precisely for 1 < z < 3. Therefore, both changes of sign occur at 
the same time and there is no sign-flip at all. The final baryon asymmetry using the full 
Boltzmann equations is in this case two (four) times bigger than our (the effective e(T)) 
approximation. For these values of K, the thermal motion of the sneutrino is important, 
so our analytic approximation gives a better agreement, specially at high T. 

Finally, for K = 10, we find no sign-fiip of the baryon asymmetry in the full Boltz- 
mann equation calculation, nor in the approximated ones, since the right-sign asymmetry 
generated in the sneutrino decay (which occurs at earlier times) compensates the sign-fiip 
due to the opening of the bosonic decay channel. As we anticipated, in this strong washout 
scenario both approximations are very good, since, in this regime, kinetic equilibrium is 
rapidly reached (Fig. V, third row), and the thermal motion of the sneutrino is negligible 
for the relevant values of z, so the effective e(T) approximation gives also accurate results. 

4. Conclusions 

The integrated Boltzmann equations are generally used to estimate the produced baryon 
asymmetry through the leptogenesis mechanism; however, these equations are not always 
a good approximation. In [3], it was shown that in the strong wash-out regime the results 
may differ by 15 — 30%. In our work, we have extended this study to two different scenarios 
where the use of the full Boltzmann equations may be relevant: the washout of a preex- 
isting asymmetry (for instance generated by a heavier right handed neutrino, and soft 
leptogenesis. 
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We consider a simplified picture which includes only the heavy right-handed neutrino 
{Ni) decays, inverse decays, and resonant scattering. This approximation is appropriate for 
T < 10^^ GeV, when the off-shell 2 2 scattering mediated by has only small effects. 

A lepton asymmetry produced during the next-to- lightest right-handed neutrino {N2) 
decay, in general is washed-out later by the lightest singlet neutrino (iVi) interactions. Then 
one wonders if a sufficient asymmetry can survive to explain the observed baryon asym- 
metry. Our results show that the wash-out obtained using the full Boltzmann equations 
is stronger than the one obtained in the usual Maxwell-Boltzmann and kinetic equilib- 
rium approximation. This result remains valid when flavour effects and different flavour 
asymmetries are taken into account. 

The second scenario that we have explored is soft leptogenesis, where the lepton asym- 
metry is produced by sncutrino interactions, and the source of lepton number and CP 
violation arc the soft SUSY breaking parameters. The lepton asymmetry produced within 
this mechanism is a pure thermal effect, because at T = the asymmetry in leptons cancels 
the one in sleptons, and only when T 7^ 0, the different statistics of bosons and fermions 
lead to a non-zero asymmetry. In this context, the use of the full Boltzmann equations 
is mandatory, since in the usual Maxwell-Boltzmann approximation the produced lepton 
asymmetry is exactly zero. However, there are various simplified approaches available in 
the literature, always assuming kinetic equilibrium for the heavy sneutrino. A possibility 
is to use an effective, temperature dependent CP asymmetry, e(r), which captures the 
main thermal effects but neglects the thermal motion of the heavy sneutrino [6]. An- 
other approach is to expand the Bose-Einstein and Fermi-Dirac distributions about the 
Maxwell-Boltzmann one, keeping only the leading order terms [9]. We have found an im- 
proved approximation, Eq. (3.49), which describes accurately the temperature dependence 
of the CP asymmetry, including the thermal motion of the sneutrino. We have compared 
our results from the full Boltzmann equations with these approximate solutions. 

We find that, in the strong washout regime {K ^1), the results obtained using the 
full set of Boltzmann equations and those obtained with the approximate equations are 
in very good agreement, because in both cases the heavy-species distribution is very close 
to equilibrium in the relevant range of z. However, in the weak washout regime (K < 1), 
there are important effects. Using the full kinetic equations, wc find that the heavy-species 
low-momentum modes are more efficiently populated, leading to a lepton asymmetry that 
can be either enhanced (as in the case K = 1) 01 suppressed (K = 0.1) with respect to the 
values obtained using the approximate equations. In fact, the difference bewteen the full 
kinetic equations and the different approximations can be up to one order of magnitude. 

Note added 

After the present paper was submitted, a new analysis of standard leptogensis with the full 
Boltzmann equations appeared [31]. The authors of this work extend previous studies by 
including, in addition to decays and inverse decays, scattering processes of the right-handed 
neutrino with the top quark. 
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A. Full BEs in standard leptogenesis 

In this appendix, we describe in detail the derivation of the full Boltzmann equations 
(BEs) relevant for standard leptogenesis, including the subtraction of the on-shell two- 
body scattering amplitude which leads to the required Sakharov condition that no lepton 
asymmetry can be generated in thermal equilibrium. 

The Boltzmann equation for the phase-space distribution of the particle species a can 
be written as: 

where the collision integral is given by: 

C[/»]= E [ dnxdnYi27r)H\pa+Px-PY) 

X [/„/x(l±/y)|A(aX^y)|2-/y(l±/„)(l±/^)|A(y^aX)|2] . (A.2) 

The sum runs over all allowed processes aX ^ y, where X and Y are multiparticle states, 
and we have used the abbreviations: 

fx,Y= n (i±/x,y)= n (A-4) 

6gX,Y beX,Y 

where -|- is for bosons and — for fermions. 
Then, for the heavy neutrino we have: 

^ - PN^H^ = dPLdpH{2n)H(^^ (pN, -PL- Ph) (A.5) 

X {fHfL{l - fN,) \A{LH N,)f - /jv,(l - + fH) \A{Ni ^ LH)\'' 
+ fHfL{^-fNj\A{LH^Ni)\'-fN,{l-h){l + fH)\AiNi^LH)f] , 

where dpx = d^px /^Exi^irf . Using CPT invariance, we can write the neutrino decay am- 

I '12 

plitudes in terms of the CP conserving total decay amplitude and the CP asymmetry 
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< 1 

Mr AT.- r,m|2 = \a(lh ^ ma\'^ = 



\A{Ni ^ LH)f = \A{LH ^ iV,)|' = ^ |^'^|' , (A.6) 



\A{Ni ^ L^)|' = \A{LH ^ iVi)|2 = 1^ 1^^'' 



We work in the approximations of thermal equihbrium for the Higgs field and kinematic 
equilibrium for the SM leptons, as described in Sec. 2, Eqs. (2.5) and (2.6). Then, substi- 
tuting the above decay and inverse decay amplitudes and neglecting terms of order 0{€fc) 
(where fc = fL — ft) of order 0{{ii/TY) (where fi is the lepton chemical potential, 
see Eq. (2.8)), we get 

' X [f/fl\l - /^J - Ml - fl'){l + /^^)] . (A.7) 
We can perform part of the phase space integration and finally obtain: 

(A.8) 

Although we do not write the arguments of the distribution functions for simplicity, energy 
conservation implies that /^'^ = /^'^(-EAr. — Eh) in Eq. (A.8). 

Regarding the evolution of the lepton asymmetry, besides the right-handed neutrino 
decays and inverse decays we need to consider the AL = 2 scattering terms, in order to 
obtain the out-of-equilibrium condition, i.e., that no lepton asymmetry is generated in 
thermal equilibrium. So we have: 

^ - PlH^^ = a - a - 25 , (A.9) 

where: 

Di = ^J dpN,dpH{2'K)H^'^) {pN, -PL- Ph) (A.IO) 

X {/iv,(l - /l)(1 + ) \A{N, - LH f - /^Vl(1 - /iv.) \A{LH N,)\'] , 

A = 2^ y" dpN,dpH{2'K)H^^^ {pM, -PL- Ph) (A.ll) 

X {/iv.(l - /z)(l + ) \A{N, ^ LH)\' - /^Vz(l - In.) \A{LH ^ iV,)|'} , 



and. 



S = j dpHdpidpH{2Tr)^S^'^'^ {pL + PH - PL - Ph) 

X {/l/h'(1 - /l)(1 + fH')\MsubiLH ^ L^)|2- 

- /z/f (1 - /l)(1 + mMsu,{LH ^ Li^)|2} , (A.12) 
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The scattering term is defined in terms of the subtracted amplitudes, since the on-shell 
contribution is already taken into account through the decays and inverse decays in the 
Di — terms. So, for example: 



with 



\Msub{LH LH)\^ = \M{LH LH)\^ - \Mos{LH LH)\^ , (A.13) 



\Mos{LH LH)f = \A{LH N,)\' ""^^j.^'^ iM^i - LSf 



I 



and analogously: 

\M„(LH ^ LH)f . i±?li l^y^I^iizig) l^y^ . (A.15) 

We define also the off-shell part, which is CP-conserving at leading order, as, 

\M,ff{LH^LH)f = \M{LH^LH)f-\A'j,f l^i^l'- (A.16) 

The width that cutoffs the resonance is Vf^, the damping rate at finite temperature: 

rf = ^ I dpLdM'i^)'s^'\pN-PL-PH) [(1 - + ) + /IVhI 1^1.1' • (A.17) 

Using the amplitudes Eq. (A.6) and the approximations in Eqs. (2.7) and (2.8), we can 
simplify the integrand of the Di — terms as follows (to first order in and fc)- 



\^D 



4* 



2 

L{/iV,(l - + )(1 + 6,) - /l(1 - /^J(l - 6,) 

-/^,(1 - /£)(! + /^^)(1 - 6,) + /^VZ(1 - /ivj(l + 6,)} 

2 - {2 [(1 - /ivJ/^Vf + /iv,(l - /r)(l + /^^)] - /a/^^ + In.)} • (A.18) 

Notice that if we only included the 1^2 processes, a lepton asymmetry would be 
generated even in thermal equilibrium, when /^Vj = and fc = 0, because the first term 
on the right-hand side of Eq. (A.18) does not vanish [32]. As anticipated, to remedy this 
problem we have to include also the subtracted contribution of the 2 <r^ 2 scattering. In 
fact, it is enough to include only the on-shell contribution to the LH LH, as the off-shell 
terms will contribute only at higher order in the Yukawas and the asymmetry. The on-shell 
scattering is [13]: 

SosiLH ^LH) = ^I dp„dpidpHi'^n)''5('\pL+PH-PL-PH) (A.19) 
X /l/^^(1 - /z)(l + ) \A{LH Ni)f ^^^^^ \A{Ni LSf , 
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Since the resonant part of the scattering term is Sqs = Sos{LH LH) — Sos{LH LH), 
substituting the decay and inverse decay ampUtudes from Eq. (A. 6) it is easy to see that it 
is order e^. Thus, within our Unear approximation we can consistently use the equihbrium 
distributions also for the leptons when calculating this term. We then rewrite the product 
of equilibrium densities as: 

(1 - fl^)il + f^^) = f%y-J^ [(1 - fl^){l + ) + fl^a] , (A.20) 

and use the identity: 

1 = J d^pN.S^^^ {PNi -PL- Ph) , (A.21) 
to obtain Sos{LH — ^ LH) at the required order in ef. 



X 

J (27r 

X j dpidpR{2^)H^'\pM^-pj^-pB)\A^^\'[{l-fl')^ . 
Notice that the integrals over the final-state particles reprodiicc^ llic thermal width of 



Ni, given by Eq. (A.17) [13]. Since / ^^2'k5{p\. - Mf) = j (2^^, we can rewrite 
Eq. (A.22) as: 

SosiLH ^LH) = ^I dpHdm.{2nr6('HpN-PL-PH)fZe^^^^^ (^^) 

(A.23) 

To this order, the Sos{LH — > LH) term is the same, just changing the (1 — eif' by (1 + e^)^, 
so altogether we find that: 

j dpHdm,{2'K)H^^\pN,-PL-PH)f'i^^e^^^'''fTr^ \^d\^ ■ (A.24) 



2El 

Using the equilibrium relation: 

we can write the Boltzmann equation for the lepton number asymmetry as: 

dfc rrdfc ^ , 1 



- PlH^ = Di-Di + 2Sos = 7^ [ dpHdpN, {2i,)^6^^) [p^^ - Pl - Ph) (A.26) 

OPL ^-C/L J 

{^^ (fN. - f2) [(1 - + ) - fi'm - IfdfH' + fNd] 
1^ r dE^. un. - a) [(1 - mil + ) - /f ] 

-Ifcifn' + fN^] > (A-27) 
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where we have integrated over part of the phase space and we have used that the heavy- 
neutrino width at zero temperature is given by Eq. (2.2). This equation does not generate 
any lepton asymmetry in thermal equiUbrium, since the right-hand side of the equation 
expHcitly vanishes in this case. Notice that in this final equation for the evolution of 
the lepton-asymmetry distribution, (A.26), we have only kept the on-shell contribution to 
the AL = 2 scattering mediated by Ni, because the off-shell part is higher order in the 
heavy-neutrino Yukawa couplings, and therefore subdominant unless these are large [33]. 

In the presence of different active lepton flavours in the plasma at T = Mj, the Boltz- 
mann equations change. The main difference is that now the right-handed neutrinos decay 
differently to the different active flavours and generate different phase space distributions 
and different asymmetries in them. Taking this into account, it is easy to see that the 
equation for the abundance of the heavy neutrino, Eq. (A. 5), now becomes: 

^ - PN.H^ = ^J dpLdpHilT^fS^''^ {PN, -PL- Ph) (A.28) 
X [fHfLA-^ - In,) \A{LaH -> iVOp - /iv,(l - /lJ(1 + fn) \A{N, ^ L^H)f 

a 

where 

\A{N, -> L^H)\' = \A{L^H ^ Ni)f = '^AaP + I^^dI ^ ^^ ^O) 

|2 _ Mi 1 2 

being yl^ ^ the CP conserving decay amplitude to the a flavour and ef the flavoured CP 
asymmetries. In the same approximations as before, i.e. up to terms of order 0{efc), we 
obtain for the heavy neutrino abundance the same equation as in the unfiavoured case, 
Eq. (A.7). 

In the evolution equations for the lepton asymmetry, we have to consider asymmetries 
for the different active flavours. The equations are now: 

- PlH^^ = Df - Df - + 5" , (A.30) 

where Df, Df are identical to Eqs. (A. 10) and (A. 11), taking into account only the 
relevant flavour in the distribution functions and decay amplitudes, i.e. — > /x,^, 
\A{N, ^ LH)\' ^ \A{Ni ^ L^H)\\ \A{N, ^ LH)\^ ^ \A{N, ^ L„i7)|', etc. 

The scattering terms, in this case, are slightly complicated as the initial and final 
leptons can have now different flavours and we have additional scatterings changing flavour 
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but not lepton number. 



(A.31) 



and 



^^ = ^1 dpHdpLdpii{2n)H^'^\pL+PH-PL-PH) X + (A.32) 



Using the subtracted amplitudes as in the unflavoured case, the on-shell part of the sub- 
tracted amphtudes, Eq. (A. 19), now is, 

Sos{LaH ^ LpH) = ^J dpHdpidpH{'^^)H^''\pL +Ph-Pl-Ph) (A.33) 

,2 TTdis - 



x/L„/f (1 - /iJ(l + ) \A{L^H 



\A{Ni^L^H)Y . 



Continuing analogously to the unflavoured case, we obtain the Boltzmann equations for 
the flavour asymmetries at first order in e and f^: 



dt 



PlH^^ = J dpHdpNi{2Tr)'^S^'^\pNi - PL - Ph) 

|2 



(A.34) 



+ 



^ j dpHdpidpH{27T)^S^^^ {pL +PH-PL- PH)fHi^ + / 



X Y{fCi3\Ma)3\lff - fCp\Map\lff - fcjMapllff - fcjMafsllff} , 

where jM^^l^jj is defined analogously to Eq. (A. 16). We have used CPT invariance, which 
implies that \M{LjiH LaH)\'^ = \M{LaH LjiH)\'^, etc., to write all the amplitudes 
with the flavour a as initial state, and also that the off-shell amplitudes are CP conserving 
at leading order in the couplings, so that |M^^|^^j = \Map\'^^^. Finally, we have summed 
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over all flavours f3 in the last line of Eq. (A. 34), since the extra contributions for (3 = a 
cancel between the first and the third terms. 

Notice that in Eq. (A. 34), we have kept also the off-shell amplitudes, even taking 
into account that they are higher order in the Yukawas. In fact, we expect |M^^|^jj to 
be of order Y^-Y^- that is indeed smaller than Y^- unless Ygj ~ 1. However notice that 
some of these substracted amplitudes contribute proportionally to fc^ and contributions 
lYaiYpil"^ fCg could be important with respect to the terms \Yai\'^fca if fca *C fcp and/or 
\Yai\ ^ \Yf3i\- Nevertheless, it is not clear a priori when these effects may play a relevant 
role. 



B. Full BEs in soft leptogenesis 



Here, we describe in detail the derivation of the full BEs relevant for soft leptogenesis. 

Using CPT invar iance and the definitions for the CP asymmetries, the decay and 
inverse decay sneutrino amplitudes can be written as: 



\A{Ni^LH)\'^ = \A{L^H^ 

|i(iVi ^ Ltilt)|2 = \A(LH^Ni)\'' 

\A{Ni^Lh)\^ = \A{Lh^Ni)\^ 

\A{Ni^Lh)\^ = \A{Lh^Ni)\^ 



, 2 

!£i| 4^12 



(B.l) 



where is the CP conserving tree-level sneutrino decay amplitude to scalars 

(fermions). 

The evolution equation for the sneutrino distribution with the same approximations 
as in Appendix A (i.e., to first order in e/. , e^. , /Li/, /Li<j) is given by 



'''' 



dt 



2% 



dpLdpH{2Tr)'^S^{pj;^, -pL- Ph) 



+ 



+ \At\' [/^Vf (1 + - + )(1 + )] } 
For the lepton and slepton asymmetries we have: 



(B.2) 



dt dpL 



- P^^^ = E (A - ^0 - 25 - ^LZt + ^^Zt - S^z + S^z + Sg , (B.3) 



dt 



P^H^^ = E (^^ - ^1) - 25 - 5^Lt + S^Li + - ~S^l + S, , (B.4) 



where 



A = j dpNidphi2TrfS^^\pNi -PL-Ph) 

X {/^^(l - /i)(l - /f )|i(iV. ^ Lhf - A(l + f^)\A{Lh - , 



(B.5) 
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Di = ^l dpN^dpHi^irfS^^^ (pN, -PL- Ph) (B.6) 
X {/iv.(l + /Z)(l + fH')\MN^ - - /^Vz(l + f},)\MLH - iV^I^} , 

and analogous expressions for Di and d\, just changing particles by antiparticles. 
As in the previous section, the scattering terms 

S = 2^ j dphdpidph{2Tr)^S^'^\pL + Ph - PL - Ph) 

- - /l)(1 - fl'mULh - , (B.7) 

= 1^ j dpHdpidpH{2Tr)^S^'^\pL + PH - PL - Ph) 

X {/Z/^'(1 + /Zt)(l + mMsubiLH ^ Zt^t)|2_ 

- /Zt/^'(1 + /Z)(l + )|M.„,(Lti7t ^ Z^)|2| ^ (B.8) 

etc., are defined in terms of the subtracted amplitudes (see eqs. (A.13),(A.14)). The thermal 
width of sneutrinos is F''^'' = rf" + T^'*'', with 

^/*' =^J dPLdPH{^^)'s^'^iPN, -PL -ph) [(1 - - /D + /r/r] > (b.9) 

rl''' = ^J dpLdpH{2^f5^'\p^-pL-pH) [(1 + /f )(1 + /^^) + r^fjf\ \At\^ . (B.IO) 

Finally, the scattering terms Sg, Sg correspond to the fast MSSM gaugino interactions 
LL •t-^ LL and are given by 

Sg = j dpidp2dp3{2Tr)'^d^'^\pL + Pi - P2 - Ps) 

X {/z(P2)/z(P3)(l - - fL{pi))\M{LL ^ LL)f- 

- fL{PL))fL{pi){l + /z(P2))(l + fz{p3))\M{LL ^ ZL)|2} 

— {particles — antiparticles} , (B-H) 

and analogously for Sg. Since these interactions are in equilibrium, we will not include 
them in our set of Boltzmann equations, but we shall impose that the chemical potentials 
for leptons and sleptons are equal, fig = fif. 

The derivation of the out-of-equilibrium condition is somehow lengthy, but completely 
analogous to the standard case described in appendix A, so we do not give many details here 
(see for example [10]). The basic point is that at 0{e), we can approximate the (s)lepton and 
anti-(s)lcpton distributions by the equilibrium ones, and then use the following relations 
between them 

(1 - 4^^)(1 - /f ) = e^-ifJ [(1 - 4^^)(1 - /f ) + , (B.12) 

(1 + )(1 + fl^) = e^-^fl \{l + f%'){l + /f") + f^Vf] , (B.13) 



L ' ' •'■« •'L 
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to reproduce the sneutrino thermal width, following the same procedure as in appendix A. 
Finally, we obtain: 

dfc 



dt 



X f/,(/iv,-/|)[(i-/f)(i-/r)+/r/i 



(B.14) 



dt 



PlH^ =^^2^1 dpj^dpHi27r)^Sip^^ -Pl-Ph)\A. 

i 

X {^s.if^, - /|) [(1 + /F)(l + fH') + fH'fl 



S|2 
i I 



C. Analytic approximation in kinetic equilibrium 



(B.15) 
(B.16) 



Our starting point are the Boltzmann equations for the (s)lepton asymmetries (3.47) and 
(3.48), with the constraint Hf = fig = fJ- imposed by fast gaugino interactions. In this 
approximation, we keep the phase space and statistical factors, crucial in soft leptogenesis, 
but we approximate the Fermi and Bose distributions by the Maxwell-Boltzmann one, since 
this is enough to obtain a non- vanishing CP asymmetry and allows to perform analitically 
the energy integrals. After integrating over the (s)lepton energy, we get 



dYc 
dz 



27r2s 



dE^e 



-En 



. N 



1 



A(l,XL,x,,)yiv(l + 2e-^^) 



dz 



X{l,XL,Xh)yNl^T'^} 

2^ 



(C.l) 



dENe-^^ le 



' N 



X{l,xj,XH)yN{l + 2e-^^) 



+ ^(,-ENil-xz+XH)/2 ^ ^-Em{1+x^-xh)/2^ ^^yMHi,x^,XH)/2 _ ^-yMHi,x^,XH)/2^ 

- X{l,X^,XH)yNtJ'T^} (C.2) 

The washout term in the above equation is the standard one, for a non-vanishing (s)lepton 
mass, so, after integrating over En gives the Bessel function /Ci (z). Notice that the chemical 
potential fi is related to the total lepton asymmetry Y^^, = + by 

r-2 



IT 



2 milC2{mL/T)+mjrlC2imi/T) 



Ycr- 



Using that 



^ ^ ^ n=0 



(yjvA/2) 



2n+l 



(2n + l)! ' 



(C.3) 



(C.4) 
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we can perform the integral over Ei\f, order by order in yjv- The resulting integrated 
Boltzmann equation for Y/^^ can be written as: 

with Fi{z) = Ff{z) + F^{z) and 

Ff{z) = X{1,x^,xh) {Xil,xj^,XH)[)Ciiz) + ICi{2z)]z 

+f-![z{3 -x^ + xh)/2] + /f [2(3 + x^- xh)/2]} , (C.6) 

F{{z) = {1 - XL + Xh)X{l,XL,Xh) {-X{l,XL,Xh)[}Ci{z) + )Ci{2z)]z 

+f([z{3 -XL + Xh)/2] + //[z(3 + XL- Xh)/2]] . (C.7) 
The function /f is given by 

I / y \ Tl-rJ. 

/I'M = E ^ (-) [A(1,X£,xh)]^"+^/C„+iM 

n=0 

Z 1 /Z\ 2 

~ -A(l, zr, XH)/Ci(2:a) + o ( - ) a:;r, a:;//)]^/C2(-2;a) 

+^Q'[Mi,^z,^/^)]'^3M + ... , (C.8) 

and // has the same structure, changing L ^ L and H ^ h. 

When thermal masses are neglected, Ff = f[, and the two different arguments of /^^ 
became the same (3^/2), leading to a global factor of 4 and the simplified function F^^\z) 
of Eq. (3.50). 

Finally, the function F2{z) is defined as 

p . N _ 97.2 [A(l, XL,Xh)f{l -XL- Xh) + [A(l, xi.xh)? 

mllC2{mL/T)+mllC2{m~jT) ' ^^'''^ 
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